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Cn ' ABSTRACT 

\o ■ 

T-H Aims. Finding tracers of the innermost regions of prestellar cores is important for understanding their chemical and dynamical evolution before 
the onset of gravitational collapse. While classical molecular tracers, such as CO and CS, have been shown to be strongly depleted in cold, 
, dense gas by condensation on grain mantles, it has been a subject of discussion to what extent nitrogen-bearing species, such as ammonia, are 
affected by this process. As deuterium fractionation is efficient in cold, dense gas, deuterated species are excellent tracers of prestellar cores. A 
comparison of the spatial distribution of neutral and ionized deuterated species with the dust continuum emission can thus provide important 
f~ i insights into the physical and chemical structure of such regions. 

Methods. We study the spatial distribution of the ground-state 335.5 GHz line of ND2H in LDN1689N, using APEX, and compare it with the 
Q distribution of the DCO^(3-2) line, as well as the 350 yum dust continuum emission observed with the SHARC 11 bolometer camera at CSO. 
» I Results. While the distribution of the ND2H emission in LDN1689N is generally similar to that of the 350 yum dust continuum emis- 
sion, the peak of the ND2H emission is offset by ~10" to the East from the dust continuum and DCO* emission peak. ND2H and ND3 
: share the same spatial distribution. The observed offset between the ND2H and DCO^ emission is consistent with the hypothesis that 
IL" ' the deuterium peak in LDN1689N is an interaction region between the outflow shock from IRAS 16293-2422 and the dense ambient 
. 5^ gas. We detect the J = 4 —> 3 line of H'^CO* at 346.998 GHz in the image side band serendipitously. This line shows the same spatial 
. distribution as DCO^(3-2), and peaks close to the 350 yum emission maximum which provides further support for the shock interaction scenario. 

_ _ I 
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1. Introduction 

Deuteration of nitrogen compounds, such as ammonia and 
N2H^, is spectacular in a number of environments including 
dark clouds, such as L134N (Tine et al. 2000; Roueff et al. 
2000; Roueff et al. 2005), low-mass star forming regions and 
prestellar cores such as LDN1689N and Barnard 1 (Loinard et 
al. 2001; Gerin et al. 2001; Lis et al. 2002a). The discussion of 
the respective contributions of grain and gas-phase processes 
in the deuteration is active, but no definite solution is yet avail- 
able. The lack of detection of deuterated water in ices toward 
low-mass young stellar objects (YSOs) by Dartois et al. (2003) 
and Parise et al. (2003) suggests that "another mechanism than 
pure solid state chemistry may be active to produce very high 
deuterium enrichment in the gas phase." 

Further constraints on the deuterium fractionation mech- 
anisms are provided by the spatial distribution of deuterated 
species, as compared with the molecular gas distribution traced 



by the submillimeter dust continuum emission. Whereas maps 
of singly deuterated species have been published in many cores 
— the relatively large line intensities of e.g. DCO^, N2D^ al- 
low easy mapping with state of the art detectors — mapping 
multiply deuterated species has proven to be a challenge, given 
the relatively low line intensities. Ceccarelli et al. (2001) have 
shown that the D2CO emission is extended around the class 
protostar IRAS 16293-2422. Roueff et al. (2005) present a 
limited map of the ND3 ground-state transition at 309 GHz 
with the CSO in LDN1689N, which unfortunately suffers from 
rather poor pointing accuracy. Because the ground state ND2H 
lo,i -Oo,o lines, at 335.5 for the ortho species and 335.4 GHz for 
the para species (See Coudert and Roueff (2006) for NH3 and 
its isotopologues line frequencies), are relatively strong (0.6 K 
in LDN1689N; Lis et al. 2006) and the atmospheric transmis- 
sion is good at this frequency, we have carried out the first map 
of a doubly deuterated species in a dense core. 
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2. Observations 

The observations have been carried out using the APEX- 
2a 350 GHz receiver of the Atacama Pathfinder Experiment 
(APEX We combined data obtained during two periods, 
August 2005 and October 2005. The receiver was tuned in 
DSB with the ND2H fines in tfie lower sideband. Tfie back- 
end was tfie facifity MPIfR Fast Fourier Spectrometer We ob- 
served a total of 10 positions towards LDN1689N (a2ooo = 
16:32:29.470, ^2000 = -24:28:52.60). Tfie pointing was cfiecked 
regularly and found accuracte to ~ 2 - 3". Tfie observations 
fiave been taken in tfie position-switcfied mode witfi a refer- 
ence position located 240" East of tfie source for tfie August 
2005 data. Tfie October 2005 data fiave been taken witfi a ref- 
erence position offset by 10' in azimutfi. Line and continuum 
maps sfiown in Lis et al. (2002b) and Stark et al. (2004) sfiow 
tfiat tfie source emission rapidly drops towards tfie East. 

Tfie data fiave been corrected for tfie sideband gain and 
tfie APEX main beam efficiency of 0.7 at 335 GHz. Tfie tem- 
perature scale used in tfiis paper is tfie main beam brigfitness 
temperature scale. Overall, tfie data calibration agree well witfi 
tfie CSO spectrum presented by Lis et al. (2006). Tfie FWHM 
beam size of APEX is 18" at 335.5 GHz. 

3. Results 

Fig. [2 sfiows tfie average spectrum integrated over tfie map, 
wfiile contour maps are sfiown in Fig.|2 Tfie ortfio- and para- 
ND2H( 1 1 1 -Ooo) lines are clearly detected witfi a mean intensity 
ratio of 2: 1, reproducing tfie ortfio-to-para statistical weigfit ra- 
tio (Lis et al. 2006). We fiave also detected tfie H'^CO+(4-3) 
line at 346.998 GHz in tfie image side band, wfiicfi is discussed 
below. 

3.1. ND2H spatial distribution 

Fig. 12 sfiows contours of tfie main ND2H emission at 335.5 
GHz (in wfiite) overlayed over a 350 yum continuum image. 
Tfie ND2H emission is clearly extended and sfiows an elon- 
gated sfiape generally similar to tfie dust continuum emission. 
Tfie limited ND3 map sfiown by Roueff et al. (2005) sfiows tfie 
same elongated pattern. Over tfie limited extend of tfie ND2H 
map, tfie line and continuum emissions are positively correlated 
witfi a correlation coefficient of 0.29. Tfiis ratfier weak correla- 
tion is due to tfie difference in peak positions of tfie ND2H line 
emission witfi respect to tfie continuum peak, tfie ND2H emis- 
sion beeing offset by ~10" to tfie East of tfie dust contimuum 
emission maximum. Tfiis offset far exceeds tfie pointing uncer- 
tainties (-2-3"). 

Tfie rigfit panel in Fig.|2sfiows an overlay of tfie ND2H con- 
tours on a DCO+(3-2) map (Lis et al. 2002b). Tfie ND2H peak 
is clearly offset from tfie DCO^ peak too. D2CO (Ceccarelli 
et al. 2000), ND3 (Roueff et al. 2005), NDjH (Loinard et 
al. 2001, Roueff et al. 2005), DCO+ (Lis et al; 2002b), afi 

' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory and the Onsala Space Observatory. 



appear at a bluesfiifted velocity compared to tfie cloud en- 
velope : Vlsr = 3.4 - 3.6 kms"' versus Vlsr = 3.8 - 4.0 
kms-^ for C'**0 and C'^O (Stark et al. 2004). Tfiis velocity 
difference was first recognized by Lis et al. (2002b) wfio sug- 
gested tfiat it results from tfie sfiock interaction of tfie power- 
ful blue lobe of tfie IRAS 16293-2422 outflow witfi tfie dense 
core. Tfie blue lobe of tfie molecular outflow is interacting and 
compressing LDN1689N, creating dense and cold post-sfiock 
gas, bluesfiifted relative to tfie cloud envelope. Tfie deuterated 
species preferably sample tfiis cold and dense material as deu- 
terium fractionation is more efficient at low temperatures and 
fiigfi gas densities. Tfie large-scale CO and continuum maps 
presented by Lis et al. (2002b) and Stark et al. (2004) present 
convincing evidence for tfiis scenario. 

Tfie spatial offset between DCO^ and ND2H provides ad- 
ditional support for tfie C-sfiock fiypotfiesis. Tfie DCO^ ion 
is detected afiead of tfie neutral species ND2H, closer to 
IRAS 16293-2422. Tfie spatial separation, 10", or -0.01 pc, is 
consistent witfi C-sfiock models witfi a pre-sfiock density of 
10"^ cm""* and a magnetic field of 30 fiG (Lesaffre et al. 2004), 
wfiicfi predict a total sfiock size of 0.04 pc. Tfie deuterium 
cfiemistry in sfiocks fias been previously studied by Pineau des 
Forets et al. (1989) and Bergin et al. (1999) in different con- 
texts. Because of tfie sensitivity of deuterium fractionation to 
botfi tfie gas temperature and tfie molecular depletion, it is ex- 
pected tfiat DCO^ will be reformed more rapidly tfian ND2H 
in tfie post-sfiock gas. Wfiile DCO^ is formed as soon as tfie 
gas temperature drops, efficient ND2H formation requires botfi 
cold temperatures and significant CO depletions, wfiicfi are ex- 
pected to occur downstream (Bergin et al. 1999). A spatial off- 
set could tfierefore exist between abundance peaks of DCO^ 
and ND2H, as we fiave detected. More detailed sfiock models, 
would be able to furtfier test tfie validity of tfiis scenario. 

3.2. ND2H abundance 

We fiave fitted tfie ND2H spectra using tfie fiyperfine fitting 
command in CLASS. Altfiougfi tfie signal-to-noise ratio is not 
fiigfi enougfi for securely deriving tfie line opacity at all posi- 
tions, it is clear tfiat tfie maximum opacity coincides witfi tfie 
maximum signal. Tfie opacity derived for tfie average of all 
positions excluding (10",0") is 0.2 + 0.48, and tfie coiTespond- 
ing excitation temperature is 8 + 3 K, in good agreement witfi 
Lis et al. (2006). Tfie spectra sfiow some indication of fiigfier 
line opacity, and lower excitation temperature (t = 3.1 1 ± 1.5; 

f = 5 ± 1 towards tfie two positions witfi fiigfiest intensities 
(0,-10") and (0,-20"). 

Tfie ND2H column densities fiave been derived assuming 
optically tfiin emission, LTE and using a uniform excitation 
temperature of 5 K. Tfiey are listed, togetfier witfi tfie H2 col- 
umn density derived from dust continuum measurements in 
Table [J We use tfie 350 jum dust continuum map obtained 
witfi SHARC II at tfie CSO (Fig. |3, and a a dust opacity of 
K350 - 0.07 cm^g"' corresponding to /cooo = 0.005 cm^g"' 
(Ossenkopf & Henning 1994) for a dust emissivity index fi -2, 
and a dust temperature of 16 K (Stark et al. 2004). Tfie 350 jim 
continuum fluxes are in fair agreement witfi tfie values reported 
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Fig. 1. Left : Full spectrum averaged over the map. p-ND2H(lii - Ooo) at 335.446 GHz, o-ND2H(l n - Ooo) at 335.514 GHz and 
Hi3CO+(4-3) at 346.998 GHz (in the image side band) are detected. Right : Enlargement showing a fit of the hyperfine structure 
of the 0-ND2H line. The temperature scale is T^b (K) and the velocity scale is relative to the LSR. 
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Fig. 2. Left: Color image of the 350 /im dust continuum emission obtained with SHARC II at the CSO, convolved to 20" angular 
resolution. The intensity scale is given in Jy/beam. White contours show the integrated intensity of the 0-ND2H line at 335.5 GHz, 
with contour levels from 0.25 to 0.45 Kkms"' with 0.05 Kkms"' spacing. Black contours show the integrated intensity of the 
Hi3CO+(4-3) hne at 346.998 GHz, with contour levels from 0.25 to 0.87 Kkms ' with 0.125 Kkms ' spacing. The map center 
is Q'2000 = 16:32:29.470, 620m = -24:28:52.60. Right: The same contours plotted on a DCO+(3-2) map (Lis et al. 2002b). 



by Stark et al. (2004), but on the lower side: we detect a max- 
imum intensity of 16 Jy in a 20" beam, to be compared with 
19.4 + 2.3 Jy in Stark et al. (2004). This dicrepancy may be 
an indication that some low-level entended emission is filtered 
out in the SHARC II image, which has been obtained in the 
"AC-biased mode", without the secondary chopper. The result- 
ing H2 column densities reported in Table[2could therefore be 
affected by this 35% uncertainty. 

The H2 column density ranges from 3.1 x 10^^ cm"^ to 
7.2x10^^ cm"^ across the map. The corresponding ND2H abun- 
dances varies from 0.7 x 10"^ to 2.6 x 10"^ relative to H2, with 
a mean value of 1.8 x 10"'^, in good agreement with previous 



work (Lis et al. 2006). ND2H appears to be remarkably abun- 
dant in this dense core. For the three positions where reliable 
ND3 data exist (Roueff^ et al. 2005), the [ND3]/[ND2H] abun- 
dance ratio is 0.01 - 0.02, and seems to increase with increasing 
ND2H column density. 

3.3. 

A relatively strong line is detected at the lower end of the ve- 
locity range in Fig.^ Careful tests have been performed using 
the CSO and APEX telescopes for identifying this line. It turns 
out that the line can be assigned to the 7 = 4 — > 3 transi- 
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Table 1. 



Position V <5V NCNDjH)' NCNDj)^ NCH,)^ 

(arcsec) kms"' kms"' 10'^^ cm"- 10'^ ctiT^ 10^^ cm"^ 

-10,10 3.63 ± 0.02 0.33 ± 0.03 9.8 ± 2 4.1 ±.5 

10,0 3.53 + 0.02 0.20 ±0.07 3.6 ±1 3.1 ±.5 

0,0 3.56 + 0.02 0.40 + 0.05 12.7 + 2 1.5 5.1 + .5 

-10,0 3.62 + 0.02 0.39 + 0.05 10.0 + 2 5.9 ± .5 

-20,0 3.75 + 0.1 0.8 + 0.3 8.9 + 2 5.4 + .5 

0,-10 3.60 + 0.02 0.30 + 0.05 15.8 + 2 5.8 + .5 

-10,-10 3.70 + 0.05 0.55 + 0.1 11.7 + 2 1.8 7.1 + .5 

-20,-10 3.55 +0.02 0.26 + 0.05 4.6+ 1 5.7 + .5 

0,-20 3.69 + 0.02 0.45 + 0.06 12.7 + 2 2.6 5.2 + .5 

-10,-20 3.72 + 0.03 0.58 + 0.08 8.9 + 2 6.5 + .5 



' Derived assuming LTE and r„ = 5 K. ^ Roueff et al. (2005).' Derived assuming Tdust ~ 16 K and k^^^q — 0.07 cm^g ^. 



tion of H'''CO+ at 346.998 GHz seen in the image side band. 
At APEX, when the system is tuned to ND2H, the H'^CO+ Hne 
falls 16MHz outside the IGHz wide IF band, but is aliased back 
into the signal band. The band edge IF filter attenuates the line 
intensity, therefore the correct intensity scale has been estab- 
lished by comparing the original data with additional measure- 
ments taken with APEX on May 22nd, 2006 with the H'^CO+ 
line centered in the signal band. The 2005 data have been mul- 
tiplied by 2.5 to match the properly tuned 2006 data. We have 
checked that the spatial distributions of H'^CO^(4-3) agree rea- 
sonably well between both datasets. H'^CO^(4-3) has a more 
sharply peaked spatial distribution than ND2H (Fig.|2jl, which 
is very similar to the DCO^(3-2) map obtained at the CSO (see 
right panel of Fig.|2j. 

4. Conclusions 

1. We show that the ND2H emission is extended in the 
LDN1689N dense core. The abundance relative to H2, de- 
rived from LTE, is ~ 1.8 x 10 which makes ND2H 
a remarkably abundant molecule in this dense core. The 
[ND3]/[ND2H] abundance ratio is 0.01 - 0.02. 

2. The ND2H emission is spatially offset from the dust contin- 
uum DCO+(3-2) and Hi3CO+(4-3) emission by ~ 0.01 pc, 
and appears to be blueshifted relative to the cloud envelope. 
Both findings are qualitatively consistent with the scenario 
of the formation of the LDN1689N dense core by the inter- 
action of the blue lobe of the IRAS 16293-2422 molecular 
outflow with the ambient material. Detailed models could 
further test this scenario. 
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